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Summary. — Accurate predictions for both signal and background events at the 
LHC are of paramount importance in order to confirm even the smallest deviations 
from Standard Model predicitions. Next-to-leading order Monte Carlo event gener- 
ators are an essential tool to reach that goal. Concerning the charged Higgs boson, 
NLO calculations of the production cross section already exist. Reiterating the cal- 
culation using a subtraction formalism enables us to implement the cross section 
into Monte Carlo generators, which can then be used by experiments. 
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1. Introduction 

MC@NLO [1] is widely used by the experiments at the LHC for top production 
cross section estimations and will thus also be the default generator for charged Higgs 
boson production, implying that this implementation had to be absolute priority [2]. 
However, MC@NLO is known to have some drawbacks for experimental use. So far, 
MC@NLO is interfaced only with HERWIG, and using another parton shower (PS) is 
not straightforward. In many cases the preferred PS is Pythia, since it is easily tunable 
to data. Another complication, which arises when using MC@NLO, is the possibility 
to encounter negatively weighted events due to the particular combination of Monte 
Carlo generator and next-to-leading order (NLO) calculation. This can be an issue if the 
experimental analysis is performed via trained multivariate techniques, which are unable 
to handle events with negative weights. A remedy to both issues can be found in the use 
of POWHEG [3]. 

Providing an implementation in two independent generators is very useful. One must 
never forget that, since event generators are complicated, man-coded computer codes, 
they are never free of bugs. At least two estimations from different codes can give a 
clearer hint as to potential errors and should be performed whenever possible. 
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2. Charged Higgs boson production at next-to-leading order (NLO) at the 
LHC 

The first NLO calculations of charged Higgs boson production in association with a 
top quark have been performed almost a decade ago [4, 5] . We have used them to check 
our results. A NLO cross section is very important to evaluate production cross sections, 
but also for tests of the universality of the charged Higgs boson couplings [6] . 
The following paragraph sketches the difficulties encountered when adding up different 
contributions of the NLO cross section. Initial-state collinear divergencies are cancelled 
with counterterms coming from the NLO PDFs. The remaining divergencies cancel 
between the virtual a v and the real cross section a R , so that the total NLO cross section 
is finite. Although the sum is finite, we need to separate these pieces in order to perform 
the integration, since they involve different final state phase spaces, 

(1) a NLO = [ da R + [ d* v . 



Various solutions to this problem have been proposed. The most popular one in the 
early days of NLO calculations was the aforementioned phase-space slicing method. A 
somewhat different approach is the so-called Catani-Seymour dipole formalism [7]. The 
general philosophy of this method is that in order to have a numerically integrable cross 
section, an auxiliary term da A is defined. This terms has two special features, namely it 
exhibits the same pole structure as da R and can thus act as a local counterterm, and it 
is analytically integrable over the singular one-particle subspace. The right hand side of 
Eq. (1) can thus be rewritten as 



(2) 
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and the integrations can now be performed over finite quantities. 

More information on this and on the implementation into MC@NLO can be found in [2]. 

3. Implementation in POWHEG 

This section briefly summarizes ongoing work concerning the implemenation of the 
NLO cross section calculation into POWHEG. The following normalised plots show a 
comparison between the pure NLO calculation and the POWHEG result when coupled 
to HERWIG for a charged Higgs boson mass of m H - = 300 GeV, a coupling to fermions 

(up-type Ui and down-type dj) of type II C cx H + Ui(^^PL + m dj tan f3P R )dj with 

Pr/l = (1 ± 7 5 )/2 and tan/3 = 30, and a center-of-mass energy of the LHC of y/S = 14 
TeV. The left plot of figure 1 displays the transverse momentum distribution of the top 
quark and the charged Higgs boson. The pure NLO curve is negative for the first bin and 
then reaches very high values. This typical behaviour is seen to be smoothened by the 
PS in POWHEG. A resummed calculation would also be similar to the PS behavior. The 
second plot shows the rapidity distribution of the top quark and the charged Higgs boson 
pair. Agreement between the pure NLO calculation and the POWHEG output coupled 
to the PS can be seen. Finally, the last plot displays the azimuthal angle between the 
top quark and the charged Higgs boson. Again, the PS regularizes the behavior of the 
NLO calculation at 6 — n. 
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Fig. 1. - Comparison between the purely NLO calculation (solid line) and POWHEG+HERWIG 
(dashed line) for the pair transverse momentum \pr top + pr _ |, the pair rapidity y top +H- arm 
the azimuthal angle <f> between the top quark and the charged Higgs boson. 



4. Outlook 

The implementation in POWHEG for large charged Higgs boson masses is complete 
and the low-mass (mcH < m t) case is in progress. However, there are still some issues 
that could be explored further. 

Our calculation relies on specific simplifications, as for example, neglecting the bottom 
quark mass with respect to other kinematic variables. Keeping the bottom quark mass- 
less is referred to as the five-flavor scheme. A compraison with the four-flavor scheme 
calculation [8] would be of great interest. An analogous comparison has already been 
performed for single top production in the i-channel [9]. Also, a proper resummation 
could be done for this process in order to compare the 1ow-pt regions with the Monte 
Carlo output. 
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